##Replication code for "drivers of resemblance: explaining the conversion of pre-electoral coalitions
##into coalition cabinets"
library(QCA)
library(tidyverse)
library(ragg)
library(SetMethods)
library(showtext)

rob_sec_sec <- read.csv("dpr.csv", sep=";")  %>%
               filter(type != "electoral") %>%
               select(-type, -legislative_type, -inter_party_primaries,
                      -low_chamber, -upp_chamber, -status, -ce, -country)

rob_sec1 <- rob_sec_sec %>%
            column_to_rownames(var = "cases")

rob_sec2 <- rob_sec1 %>%
            mutate(lwpol = calibrate(lwpol, type = "fuzzy", method = "direct", threshold = "e=3, c=2, i=1")) %>%
            mutate(hlpol = calibrate(hlpol, type = "fuzzy", method = "direct", threshold = "e=2, c=3, i=4")) %>%
            mutate(htemp = calibrate(htemp, type = "fuzzy", method = "direct", threshold = "e=92, c=61, i=55")) %>%
            mutate(lplp = calibrate(lplp, type = "fuzzy", method = "direct", threshold = "e=0.6, c=0.399, i=0.3")) %>%
            rename(LWPOL = lwpol,
                   MAJ = maj, 
                   HLPOL = hlpol,
                   HTEMP = htemp,
                   LPP = lplp,
                   CR = cr)
##Analysis of Sufficiency
TTsec1 <- truthTable(rob_sec2, outcome = "CR", conditions = c("MAJ", "LWPOL", "HLPOL", "HTEMP", "LPP"),
                     complete = TRUE, incl.cut = 0.8,
                     sort.by = "incl, n+",
                     show.cases = TRUE)

ESA <- findRows("~MAJ*~HTEMP", obj = TTsec1)

##The code below produces Table A.9 in the Supplementary Material. Pathways have been rearranged for better readability. 
intermediate_result <- minimize(TTsec1, include = "?", details = TRUE, dir.exp = "1,1,1,1,1", exclude = ESA)

stargazerSol(intermediate_result, 
             outcome = "CR",
             sol = 2,
             show.cases = T,
             type = "latex")
